The strongest bounds on active-sterile neutrino mixing after Planck data* 
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Light sterile neutrinos can be excited by oscillations with active neutrinos in the early universe 
and contribute as extra-radiation, parameterized in terms of the effective number of neutrino species 
N e g. This parameter has been measured to quite a good precision by the Planck satellite experiment, 
yielding N e g = 3.30 ± 0.27 at 68 % C.L. We use this result to update the bounds on the parameter 
space of (3+1) sterile neutrino scenarios, with an active-sterile neutrino mass squared splitting in 
the range (10 -5 — 10 2 ) eV 2 , in both normal and inverted mass hierarchies for the active and sterile 
states. For the first time we take into account the possibility of two non-vanishing active-sterile 
mixing angles. We find that the bounds are stronger than those obtained in laboratory experiments. 
In fact, we get active-sterile mixing angles sin 2 6n < 10~ 2 ' 5 for mass splittings Am^ > I0 _1 eV 2 . 
This result leads to a strong tension with the short-baseline hints of light sterile neutrinos. In 
order to relieve this disagreement, modifications of the standard cosmological scenario, e.g. large 
primordial neutrino asymmetries, are required. 

PACS numbers: 14.60.St, 14.60.Pq, 98.80.-k 26.35.+C 



Introduction. — Light sterile neutrinos have been ad- 
vocated as a possible solution to some puzzling results 
found in neutrino oscillations, see [lj for a recent re- 
view. In particular, m ~ 0(1) eV sterile neutrinos mixing 
with the active states have been proposed to solve differ- 
ent anomalies observed in short-baseline neutrino experi- 
ments, in the -+ v e oscillations in LSND Q and Mini- 
Boone [1,0] experiments, and in77 e and v e disappearance 
revealed by the Reactor Anomaly p| and the Gallium 
Anomaly [g, respectively. Scenarios with one (dubbed 
"3+1") or two ("3+2") sub-eV sterile neutrinos 0, i] 
have been proposed to fit the different data. On the other 
hand, lighter sterile neutrinos with Am 2 ~ 10 -5 eV 2 can 
explain the absence of the upturn in the solar neutrino 
energy spectrum Q . 

The main theoretical motivations for these states are 
perhaps not so strong, though light sterile neutrinos with 
a sizable mixing emerge in several models (see, e.g., fioj 
and references therein). In any case, their role as mes- 
sengers of new physics is relevant enough to justify an 
open mind attitude and a close looking for any, yet tiny, 
evidence for new effects beyond the too much successful 
Standard Model. 

The hunt for sterile neutrinos in laboratory exper- 
iments is currently open. Different techniques have 
been proposed to search for these elusive particles (see, 
e -g-j lull)- Indeed, since every experimental mcasure- 
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ment has its own systematic uncertainties and its own 
recognized or un-rccognized loop holes, in order to cor- 
ner the sterile neutrino parameter space it is worth using 
as many handles as possible (see, e.g., fl2h. Cosmology 
is one of the most powerful tool (see, e.g., [l3(). Adding 
exotic contribution to the radiation content in the uni- 
verse, usually expressed in terms of the effective number 
of excited neutrino species, N e g, has in fact, a big im- 
pact on both the Cosmic Microwave Background (CMB) 
anisotropy map fl4T4l7l |. and the Big B ang Nucleosyn- 
thesis (BBN) nuclear species yields (la , [l~9| . The stan- 
dard expectation for this parameter is N e g = 3.046 [20l ]. 
If low-mass sterile neutrinos exist and mix with active 
flavors, they can be thermally excited by the interplay 
of oscillations and collisions, producing a larger value of 
iVeff. Cosmological constraints on sterile neutrinos based 
on their contribution to the extra-radiation have been 
presented in several papers, see e.g. fl3l. [2lT[23j . 

In the last few years, a possible cosmological hint of 
light sterile neutrinos (see e.g. [U, [25[) was found by 
combining the result in the best fit of WMAP, SDSS 
II-Baryon Acoustic Oscillations and Hubble Space Tele- 
scope data, yielding a 68 % C.L. range on Neg = 
4.34+0 88 Gj for a ACDM universe. The recent results 
of WMAP-9 [ll|, SPT and ACT [H, exploiting the 
damping tail features at high multipoles, have weakened 
this evidence to less than 2-cr. 

A recent breakthrough in constraining the dark ra- 
diation content in the early universe is represente d by 
the first data release of the Planck collaboration [26| . 
a satellite experiment with unprecedented sensitivity in 
the high multipole range. Indeed, one of the main result 
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of this new CMB anisotropy map is the quite accurate 
estimate of the rclativistic degrees of freedom at recombi- 
nation epoch, N e g = 3.30±0.27 at 68 % C.L., a result ob- 
tained combining Planck, WMAP, Baryon Acoustic Os- 
cillation and high multipole CMB data [26j]. Within 1-er 
this is compatible with the standard expectation, but still 
leaves room for almost an extra neutrino species at 95 % 
C.L. 

This result is the motivation of this paper, where we 
present an update of the cosmological bounds on light 
sterile neutrinos. We focus on (3+1) scenarios, con- 
sidering a broad range for active-sterile neutrino mass 
splitting, which cover the regions where laboratory hints 
emerge. We have chosen to consider a minimal scenario 
(one extra sterile state only), and not to cover the (3+2) 
case, where it seems harder to fit the neutrino mass 
bound from large scale structure [lil . This case also ap- 
pears to be disfavored by Planck results unless one con- 
siders large neutrino asymmetries to suppress the sterile 
neutrino production, see e.g. p2l. \2l\. 

Setup of the (3+1) flavor evolution. — We consider 
one single light sterile neutrino v s , which mixes with 
the active neutrino states v e ,v^,i> T . The flavor eigen- 
states v a are related to the mass eigenstates vi via a 
unitary matrix U = U(6 12 , 9 13 , 2 3, U , #24, M ) 0, Hp 
where we order the flavor eigenstates in such a way that if 
all mixing angles are vanishing we have (y ei v^, v Tl v s ) = 
( v\, V2, v%, va). In the following we fix the values of the 
three active mixing angles to the current best-fit from 
global analysis of the different active neutrino oscilla- 
tion data sin 2 6»i2 = 0.307, sin 2 6» 23 = 0.398, and 
sin 2 6*i3 = 0.0245. Concerning the mixing angles between 
active and sterile neutrinos we choose as representative 
range 10~ 5 < sin 2 9 tA < lO" 1 (i=l,2,3). 

The Av mass spectrum is parameterized as M 2 = 



diag (m 2 , mf 



Am^ , m{ + A??i 



The 



solar and the atmospheric mass-square differences are 
given by Am^j = to 2 — rr 
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| - m{ = 7.54 x 10~ 5 eV 2 and 
2.43 x 10" 3 eV 2 , respectively [H]. 



Depending on the sign of Am 2 ; and Am^ we define an 
active normal mass hierarchy (NH, 7773 > mi) or an ac- 
tive inverted mass hierarchy (IH, 7713 < mi) and a sterile 
normal mass hierarchy (SNH, 777,4 > mi) or a sterile in- 
verted mass hierarchy (SIH, 777,4 < mi)- In our study we 
consider the following ranges 10~ 5 < Am^/eV 2 < 10 2 
in SNH and 10" 5 < \Am 2 41 \/eV 2 < 10~ 2 in SIH. Note 
that in SIH larger values of lAm^l are disfavored due to 
the cosmological bounds on the neutrino masses 

nam 

The neutrino (antincutrino) ensemble in a medium, as 
in the early universe plasma, is described in terms of 
a 4 x 4 momentum-dependent density matrix g p (g p ). 
To solve the full set of momentum dependent equa- 
tions of motion [30l | turns out to be a computation- 
ally demanding task (see, e.g., [2^, I27} for recent stud- 
ies). Since our aim is to perform an extensive scan of 
the sterile neutrino parameter space, in order to carry 
out a more treatable numerical analysis, we will con- 
sider the averaged-momcntum approximation, based on 



the ansatz, g p (T) -> f FD {p)p{T) (see [3l|) where p{T) 
is the density matrix for the mean thermal momentum 
(p) = 3.15 T, and /fd{p) is the Fermi-Dirac neutrino 
equilibrium distribution, and similarly for antineutrinos. 

The evolution equation for the momentum-averaged 
density matrix p, describing the neutrino system, is the 
following 



latrix p, 

Mm 



C[p,p] 



(i) 



and a similar expression holds for the antineutrino ma- 
trix p. The evolution in terms of the comoving observer 
proper time t can be easily recast in function of the tem- 
perature T (see [3l| for a detailed treatment). The first 
term on the right-hand side of Eq. ([T|) describes the flavor 
oscillations Hamiltonian, 
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where M 2 = M. 2 U is the neutrino mass matrix, while 
the terms proportional to the Fermi constant Gf encode 
the matter effects in the neutrino oscillations. In par- 
ticular, the term is related to the energy density of 
pairs, E v to the energy density of v and D, and is 
the v — v interaction term proportional to the neutrino 
asymmetry. In the following, we will consider the most 
conservative scenario, with zero neutrino asymmetries, 
or as small as the baryon asymmetry, t\b ~ 6 • 10 - . 
Finally, the last term in the righ-hand side of Eq. (|T|) is 
the order G 2 F collisional term. 

The matter terms can induce Mikheyeev-Smirnov- 
Wolfenstein (MSW)-like resonances [32[ between the ac- 
tive and sterile states when they become ©(Am 2 ^). The 
possible resonances depend on the different mass hierar- 
chies in the active and sterile neutrino sector. In particu- 
lar, in the absence of neutrino asymmetries the resonance 
condition can be satisfied (in both neutrino and antincu- 
trino sectors) only for Am 2 ; < [311 ] . When more than 
one Am 2 ^ is negative, multiple resonance can occur, af- 
fecting the sterile neutrino production. 

Bounds on active-sterile v mixing after Planck. — The 
quantity we are exploiting to constrain the sterile neu- 
trino parameter space is the overall non electromagnetic 
radiation content, parametrized via 7V e ff, 



N eS = ±Tr{p + p) 



(3) 



Our bounds are given comparing this number with the 
one measured by Planck experiment, N e s < 3.84 at 95 
% C.L. [26j l. We present our exclusion plots in the planes 
(Am^jSin ^4). We checked that our results are not 
sensitive to the active neutrino mass ordering. Therefore, 
in the following we show our findings only for NH, while 
we consider both SNH and SIH. Since the results as a 
function of sin 2 #34 and sin 2 824 are very similar, we omit 
to present the sin 2 #34 case. In Figure 1 we show the 
exclusion plots at 95 % C.L. In the upper panels a), b) 
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FIG. 1: Exclusion plots for the active-sterile neutrino mixing parameter space for different scenarios (see the text for details). 



we consider SNH, while in the lower panels c), d) we show 
the SIH case. Left panels refer to the exclusion plots in 
the plane (Am^, sin 2 614) for different values of sin 2 024, 
while right panels refer to the plane (Am 2 1; sin 2 024 ) for 
different values of sin 2 6*14. In all the cases, sin 2 6*34 is 
fixed to zero. The excluded regions are those on the 
right or at the exterior of the (closed) contours. We now 
discuss the different panels in more details. 

Panel a) The most conservative bound corresponds to 
sin 2 024 = 0, where for Aro^ > I0 -2 cV 2 the exclu- 
sion contour is a straight line in this plane. Lowering 
the value of Am^ one triggers first a V4 — V3 resonance 
(when Am 2 ! < Am^) and then also a V4 — V2 resonance 
(when Am 2 ! < Ara^) which is the dominant one. These 
produce the changes of the slope in the exclusion plot. In 
particular, increasing the value of sin 2 024, the constraint 
on the parameter space becomes stronger. Large values 
of sin 2 024 would dominate the sterile neutrino produc- 
tion and this excludes regions otherwise permitted if only 
sin 2 0i4 were non-zero. In particular, the only part that 
remains open is the transition region between the effi- 
cient non-resonant production range at large Am^ and 
the one of resonant production at small Am^ . For com- 
parison, we also show the slice at sin 2 024 = 10~ 2 of the 
95 % C.L., allowed region obtained from the global anal- 



ysis of short- baseline oscillation data @, HU (filled region 
in the up right part of the plot denoted by SBL). We see 
that it seems to be completely ruled out by the cosmo- 
logical bound. We also plot the 90 % C.L. expected sen- 
sitivity of the KATRIN experiment (measuring the spec- 
trum of electrons from tritium beta decay) after 3-years 
of data taking j34j. Also this region would be already 
excluded by cosmology in absence of any significant pri- 
mordial neutrino asymmetry [22|, [2^, [53, Hl[ . Finally, we 
represent with a square in the lower part of the plot, the 
region of parameters corresponding to a light sterile neu- 
trino with Am^ ~ 10- 5 eV 2 and sin 2 M - ICT 4 - 10" 3 , 
suggested to solve the problem of the upturn of the solar 
neutrino spectrum 

Panel b) The description of the exclusion plot is anal- 
ogous to the one of panel a), with the roles of 0i4 and 
024 interchanged. In particular, the region of resonant 
sterile neutrino production is at Am 2 ! ~ I0 -3 eV 2 when 
a fi — U3 resonance is efficient. We have also drawn the 
slice at sin 2 9u = 10 -1 ' 5 of the 95 % C.L. allowed re- 
gion obtained from the global analysis of short-baseline 
oscillation data @, [H| , which is another view of the SBL 
region shown in panel a), and the exclusion curve ob- 
tained from the combined analysis of the data of and 
disappearance experiments. One realizes that also in 
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this case the region is excluded by the cosmological limit. 

Panels c) and d) Since Am^ < also a v^-vi reso- 
nance is present. This leads to an increase in the pro- 
duction of sterile neutrinos, with respect to the case of 
SNH. Therefore, the excluded regions in the parameter 
space are larger than the corresponding ones in the upper 
panels. 

Conclusions. — In this paper we have exploited the 
very recent measurement of N e g provided by the Planck 
experiment to update the cosmological bounds on (3+1) 
sterile neutrino scenarios under the assumption of van- 
ishing or very small neutrino asymmetries, of the order of 
the baryonic one. To our knowledge, for the first time it 
is shown how the constraints change if two active-sterile 
mixing angles are considered. 

We find that the sterile neutrino parameter space 
seems now to be severely constrained, and the excluded 
area covers the region accessible by current and future 
laboratory experiments. From the results of our analysis 
we conclude that there is a tension with the sterile neu- 
trino hints from short-baseline experiments. Notice that 
combining Planck findings with other data might fur- 
ther strengthen the bounds on N c ff . For example, adding 
in the analysis the primordial deuterium determination 
of Ref. HU , compared with the BBN theoretical expec- 
tation as function of baryon density and N e g, leads to 
N eS < 3.56 at 95 % C.L. H. This means that future 



2 H measurements reducing the present spread of different 
Quasar Absorption System results would lead to stronger 
bounds on sterile neutrino mixing parameters. 

In order to reconcile the laboratory signals in favor of 
extra sterile neutrino degrees of freedom with the cosmo- 
logical bounds one should introduce some extra parame- 
ters in the so far extremely succesful standard cosmolog- 
ical model, as for example, large neutrino-antineutrino 
asymmetries, L v = (n v — n p )/n 7 > 10~ 2 [13: HH, which 
might inhibit the sterile neutrino production in the early 
universe. After all, the fact that a completely satisfactory 
model of everything might not yet achieved is welcome, 
as it would continue to trigger the curiosity of physicists 
to look for what is, hopefully, beyond the corner. 
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